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1. Introduction 

We can obtain the wave function renormalization factor Z' defined in the RI'-MOM scheme 
and the mass renormalization factor Z q ■ Z m defined in the RI-MOM scheme from the staggered 
quark propagator using non-perturbative renormalization (NPR) method. We generate the stag- 
gered quark propagators using the momentum source in the Landau gauge on the MILC coarse 
lattices [[j], §]. Here, we present results of the data analysis after the projection. 

2. Mass Renormalization 

Let us consider a staggered fermion propagator. 

S f cd {x u x 2 ) = (x{(xi)X f Ax2)) , (2-1) 

where / is a flavor index, c, d are color indices. Here, x\ and x 2 represent the position coordinates 
on the lattice with x\ ,x 2 G Z 4 . The lattice spacing is a. 

In the normal Brillouin zone, we use p,q as the momentum, and, in the reduced Brillouin zone, 
we use p, q as the momentum as follows, 

p,qe( ,-] 4 , p,q£ (- — ,— ] 4 , q = q+7l A , p = p + 7I B , (2-2) 

a a 2a 2a 

n 

where A, B are hypercubic vectors whose element is or 1, and % A = —A. 

a 

The staggered quark propagator is defined as 

, f 

y cc ,(p + K A ,q + K B ) = (x{(p + 7l A )x c ,(q + n B )). (2.3) 
Using the Fourier analysis, one can show the following relationship. 

S f cc , (p + n A ,q + n B ) = 8^ (p - q) [S{p)] f mcc , (2.4) 

where 

8^(p) = (2a) 4 £ e'Py (2.5) 
.yeW 4 

Here, y represents a position coordinate of a hypercube whose lattice spacing is 2a, y G W 4 , and 
W denotes one-dimensional lattice whose spacing is 2a. By setting p = q, the quark propagator 
becomes 

S{. c , (p + 7l A ,p + K B ) = 8^ (0) [S{P)] f m cd = V [S(P)] f AB;cS . ( 2 " 6 ) 

where V is lattice volume factor. 

V = (2a) 4 £ 1 = L 3 x T (2.7) 

yeW 4 

where L (T) is lattice size in the spacial (time) direction. 
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To obtain the propagator, we have to solve the staggered Dirac equation. 

(fl s + m f )iy( = h 
f 1 



(2.8) 



where i is a gauge configuration index, {^ s + m^)i is Dirac operator for staggered fermion. h is a 
source vector, and \jff is a solution vector of the Dirac equation for a specific gauge configuration 



Sj {x u x 2 ) 



1 



{0 s + m f ) 



(2.9) 



where sj (x\ ,x 2 ) is a propagator for a specific gauge configuration i. A thermalized quark propaga- 
tor is defined as 



1 N 

S f ccJ (xi,x 2 ) = - ^£s{. ccl (xux 2 ) 



(2.10) 



where N is the number of gauge configurations. 

and s{ are quite different, sf do not conserve the momentum, because there are external 
gluons which live for a short time in a Monte Carlo evolution time. But conserve the momentum 
in the reduced Brillouin zone. The difference is shown in Figure [l[ 
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Figure 1: S{ has temporal external gluons, while all the gluons are contracted in S?. 
The solution vector y/f is 

Vi(*i)=« 4 E 4<fe( x i' x 2)^(^2) 

x 2 GZ 4 

We set the source vector hb(x 2 ,p + Kb) = e~ l ( p+7lB ^ X2 8t c - 

a: 2 gZ 4 



(2-11) 



(2.12) 
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After we calculate \jff, ac (x\,p + %b) using conjugate gradient method for each c, we can obtain the 
full matrix of y( ab (x\ , p + 7Tg): 



„4 N 

II 

i xieZ 4 



a N z i ^^vU^P+^^mL^ 



(2.13) 



The inverse bare propagator can be expressed as follows, 

i&W&S = K 1 + Z S )m f (r®T) AB + (1 + L V ) £ i sin (^)(^®T 
+ Z r m^ ^ sin(p M a)(sin(^ v a)) 3 (7 MV ® 1' 



,1B 



'AS 



•Ea -sin(p M a)(sin(p v a)) 3 (sin(p p c[)) 5 (7 A1V p(g)i; 



+I P m£ £ sin(p M a)(sin(p v a)) 3 (sin(p p a)) 5 (sin(p CT a)) 7 (7p V p (7 (g) l) AB ] rc 



which is derived from the lattice symmetry [Q]. The definition of (75 ® ^f) ab is given as 



1 



(75^^U = 4^[7]757b7f] 
1 

16, 



(75®^U = T rK-l) /1C (75»^U(-l) DS (2-14) 

The renormalization of propagator is Sg(p) = Z q sl(p), and the mass renormalization is defined 
by m« = Z m m,Q. Here, Z g is wave function renormalization factor for quark field, Z m is mass 
renormalization factor, niR is renormalized quark mass, mo is a bare quark mass and the subscript R 
denotes a renormalized quantity, the subscript denotes a bare quantity. Unless specified, we use 
the convention of m = mo in this paper. 
The RI'-MOM scheme prescription is 



to ^ p 



Z' 



Z m m + Ci 



Ore) 



48 



Tr^l®!)^ 1 ^)] 



(2.15) 
(2.16) 



where pp = sia(ap^) and p =Y*u.Pu.- So the renormalized propagator can be rewritten as 

x -1 



(I + E5) m R - 



Thus, we can write the Z and Z m as follows, 



(l+£y)Z„ 



Z ? = (l+Iy), Z„ 



(2.17) 



(2.18) 
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3. Results 



We generate staggered fermion propagators for 5 quark masses and 6 external momenta with 
0.5 < \ap\ < 0.75. am = 0.01,0.02,0.03,0.04,0.05. The fitting function suggested in Ref. % g, g, 
0] is 

y i = ^Tr[S- l (p)F i ] (3.1) 
f q (m,a,p)= Cl ('l+r 1 [log(/) + 2^] N ) +C2 log(/)+ C3 [log(p 2 )] 2 + C4^ (3.2) 



+ c 5 ^^log(/3 2 ) +c 6 (am) + a p 2 + c%{p 2 ) 2 +c 9 p 4 (3.3) 



2 f„ m \2 A2 



f n (m,a,p) = 4 + W +r 2 [log(p 2 ) + + log(l + -?—)]) (3.4) 



4 [log(p 2 )] 2 + J 5 ^+,/ 6 -^log(p 2 ) (3.5) 



^2 . . /-2\2 , . -4 ' 



+ d 7 p* + ds(py+d 9 p« \ (3.6) 

where the anomalous dimension T/ is 

„ a v 4 „ a, 16 

r 2 = - T ± -^r- (3-7) 



(47i) 3' 471 3 

Here, y,- represents a data point obtained by some projection P ; -. 

Let us consider a data analysis for Z' q with a vector projection: Py = (y M <g> \ )p/p 2 . We use the 
uncorrelated Bayesian method to fit the data to f q (X) by imposing the following prior condition: 
c\ = 1 ±0.5a s , C2-5 = 0±2a i 2 , and C6-9 = 0±2. Here, X represents m,p,a collectively Here, 
note that the prior information on c 1-5 comes from the lattice perturbation theory [0]. In order to 
investigate the fitting quality, let us define Ary as Ary =y v — c 9 p 4 . In Fig. 2(a) and [2(b)| , we present 



Ary and yy, respectively. The fitting results are given in Table As you can see in the plots, the 
fitting quality is quite good. 
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(a) Ar s 

Figure 3: Ars vs. p 2 (left), and y$ vs. am (right) 




Let us turn to the data analysis for the scalar projection: F$ = (1 <8> 1). We use the Bayesian 
method to fit the data to f m (X). The prior conditions are di = 1 ±0.5a s , c?3_6 = 0±2a 2 , and 
dj-9 = 0±2. Let us define Ars as ^ r s = ys — d 9 (am)p 4 . In Fig. ||, we show Ars and y^. As one can 
see in the plots, the fitting quality is somewhat poor with # 2 /d.o.f = 1.24(39) for the uncorrected 
Bayesian fitting. The fitting results are summarized in Table ||. 

In Fig. EL we define the y-axis variables as 



Z 'a{P) = fq(P'> m = °> Cl-9 = 0) 



Z q -Z m (p) = — f,„{p; m = 0, d\ = 0, d 7 - 9 = 0) 
am 



(3.8) 
(3.9) 



We estimate the statistical errors using the jackknife resampling method. As you can see in the 
plots, the minimum of statistical errors is located at \p\ = 2.084GeV for Z' q and at \p\ = 2.190GeV 
for Z q -Z m . Hence, we choose \p\ = 2GeV as our optimal matching scale. Our preliminary results 
are 

(3.10) 



z;(p = 2GeV) 


= 0.9810(46), Z q -Z m (p = 2GeV) 


= 1.0551(52). 




Table 1: Fitting results forZ^. 






Cl C3 C4 


C5 


0.931( 18) 


0.185(47) 0.161(28) -0.096( 18) 


0.263( 39) 




Cl c% eg 


X 2 /d.o.f 


-0.012( 13) 


0.97( 17) -1.17( 19) 0.373(27) 


0.35( 11) 


Table 2: Fitting results for Z q ■ Z m . 




Cl c 3 c 4 


C5 


0.02992( 38) 


1.143( 11) -0.205( 19) -0.117(20) 


-0.0346( 82) 




Cl c 8 c 9 


X 2 /d.o.f 


0.0414( 98) 


2.358( 82) -2.66( 16) -2.55( 39) 


1.24( 39) 
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We plan to cross-check these results against those obtained using the bilinear operators in near 
future. 
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